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summary 


This  Report  shows  how  a  direct-lift  motivator  can  be  used  to  alleviate  the 
response  to  vertical  turbulence  of  a  rigid,  combat-type  aircraft.  Optimal- 
control  and  parameter-optimisation  techniques  are  used  to  design  both  'open'  and 
'closed'  loop  control  systems.  Alternative  criteria  other  than  simply  reducing 
the  normal  acceleration  response  are  discussed.  The  effects  of  system  non- 
linearities  such  as  position  and  rate  limits  are  explored.  Both  discrete  and 
continuous  models  of  turbulence  are  used  in  the  analysis.  Some  nonlinear 
control  solutions  are  discussed. 
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1  INTRODUCTION 

This  Report  deals  with  the  problem  of  designing  a  control  system  incorpor¬ 
ating  direct-lif t-control  (DLC)  to  improve  the  ride  performance  of  a  modern 
fighter  aircraft.  It  is  one  of  a  series  dealing  with  the  implementation  of  poss¬ 
ible  performance  benefits  of  active  control,  as  applied  to  military  combat 
aircraft.  Active  control  is  defined  here  as  the  use  of  feedback  or  feedforward 
control  to  change  the  dynamic  characteristics  of  the  aircraft.  There  are  basic¬ 
ally  two  types  of  benefit:  firstly,  that  associated  with  pilot  handling 
characteristics  including  tracking,  weapon  aiming,  ’care-free'  manoeuvring  etc, 
and  secondly,  that  associated  with  fully  automatic  control  systems,  independent 
of  pilot  input,  including  ride  control,  gust  and  load  alleviation  and  flutter 
control.  'Active  control'  also  implies  the  use  of  any  available  motivator, 
including  not  only  the  conventional  elevator,  aileron,  flap  and  rudder,  but  also 
other  possible  motivators  such  as  suitably  placed  canards  or  deflected  thrust 
jets.  From  basic  principles,  the  more  motivators  available  to  the  control 
engineer  the  more  the  design  flexibility  and  range  of  control  he  has  over  the 
aircraft  dynamic  response. 

Research  into  ride  quality  of  aircraft  began  more  than  50  years  ago  and 
the  benefits  of  high  wing  loading  and  a  low  lift-curve  slope  were  established. 

The  advent  of  high  wing  loading  aircraft  with  their  inherently  improved  ride 
qualities  somewhat  negated  the  need  for  ride  smoothing  control  systems.  How¬ 
ever,  the  modern  requirement  for  multipurpose  aircraft  having  both  good 
manoeuvring  characteristics  at  high  altitude  (ideally  requiring  low  wing  loading 
and  high  lift  curve  slope)  and  good  ride  performance  at  low  altitude  has  revived 
interest  in  ride  control.  There  is  now  the  possibility  that  the  poor  ride 
qualities  near  ground  level  of  a  low  wing  loading  aircraft  may  be  significantly 
improved  through  the  use  of  active  control. 

A  synopsis  of  the  current  knowledge  of  ride  qualities  can  be  found  in 
Ref  1.  At  present,  few  generally  accepted  criteria  for  defining  good  or  accept¬ 
able  ride  qualities  are  available,  but  clearly  the  motion  of  aircraft  in  response 
to  atmospheric  turbulence  is  the  major  factor. 

A  pre-requisite  for  the  design  of  a  ride-smoothing  system  is  a  quantitative 

2 

specification  of  the  atmospheric  environment  .  Atmospheric  turbulence  is 
usually  modelled  analytically  as  a  statistical  process  and  aircraft  motion 
resulting  from  the  turbulence  excitation  is  generally  deduced  on  the  basis  of 
random  process  theory.  However,  since  the  standard  statistical  models  do  not 

completely  describe  atmospheric  turbulence  other  discrete  forms  of  input  must 

.  .  .3-5 

not  be  ignored  . 
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Only  the  longitudinal  axis  has  been  considered  in  this  Report  and  the 
control  system  has  been  implemented  using  two  motivators:  the  conventional 
elevator,  and  split  flaps  to  generate  DLC.  The  design  of  the  system  was  also 
constrained  as  follows :- 

(i)  all  modes  of  the  controlled  system  are  damped  to  at  least  0.5  relative 
damping; 

(ii)  the  elevator  and  DLC  motivator  rates  and  positions  are  limited  to 
within  practical  values; 

(iii)  a  representation  of  the  actual  hardware  is  included  in  the 
mathematical  model. 

In  order  to  obtain  an  insight  into  some  of  the  problem  areas,  the  design 

was  first  treated  using  a  simple  model  of  the  aircraft  dynamics  plus  first- 

order  filters  representing  the  actuation  systems,  including  both  elevator  and 

6  8 

DLC.  The  design  technique  involved  the  use  of  parameter  optimisation  *  ,  optimal 

7  Q  8  ft  9 

control'*  »  time  response  ,  frequency  response"*  and  root-locus  programs,  both 
digital  and  hybrid.  This  initial  design  approximated  to  pure  state  feedback 
and  enabled  a  broad  assessment  to  be  achieved.  These  initial  control  laws 
were  then  applied,  with  some  small  modifications,  to  a  more  complete  and  realis¬ 
tic  model  of  the  aircraft  plus  control  system,  including  representation  of 
actuators,  power  controls,  sensors  with  their  associated  noise  filters,  and 
various  nonlinearities. 

Two  flight  cases  have  been  considered,  one,  a  low-speed,  high  incidence 
case  and  the  other  a  high  Mach  number  case,  both  at  low  altitude.  Normal 
accelerations  at  various  positions  along  the  fuselage  have  been  calculated, 
although  the  optimisations  were  mainly  done  for  normal  accelerations  at  the 
centre  of  gravity.  For  the  closed-loop  systems  the  reduction  in  rms  normal 
acceleration  was  relatively  insensitive  to  position  along  the  fuselage.  Worth¬ 
while  theoretical  reductions  of' the  rms  of  normal  acceleration  have  been 
achieved  (>30%). 

Comparisons  are  made  between  results  obtained  using  three  different  types 

of  turbulence  model:  band-limited  white  noise,  the  conventional  Dryden  spectrum, 
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and  families  of  discrete  gusts  *  .  Since  the  gust  scale  length  of  the  Dryden 
spectrum  is  a  function  of  height,  this  parameter  has  also  been  varied.  Per¬ 
centage  reductions  in  rms  of  normal  acceleration  are  shown  to  be  proportional  to 
the  DLC  actuator  activity.  The  effects  of  design  constraints  due  to  actuator 
rate  and  position  limits  have  been  considered. 
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Since  it  has  been  suggested  that  an  appropriate  measure  of  'bumpiness* 
should  involve  the  number  of  peaks  or  exceedances  in  normal  accelerations  per 
unit  time,  a  preliminary  assessment  in  terms  of  this  criterion  has  been  included. 
It  has  been  shown  that,  compared  to  the  datum  aircraft,  for  closed-loop  ride- 
control  systems  the  number  of  zero  level  crossings,  or  total  number  of  bumps, 
will  increase,  although  the  number  of  g-exceedance  levels  will  decrease  above 
certain  amplitudes.  Thus  although  the  rms  of  normal  acceleration  may  have  been 
reduced  drastically,  the  number  of  g-exceedances  below  a  certain  level  will  have 
increased.  In  discrete  gust  terminology  this  means  that  the  ’tuned’  gust  length 
has  decreased  relative  to  that  of  the  basic  aircraft.  The  correlation  between 
variances  of  aircraft  states  derived  from  conventional  statistical  turbulence 

models  and  peak  response  amplitudes  obtained  using  discrete  'worst-gust' 
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theory  is  demonstrated.  The  total  number  of  bumps  is  shown  to  depend  on  the 
ratio  between  the  rms  of  rate  of  change  of  normal  acceleration  and  the  rms  of 
normal  acceleration  itself.  Thus  some  on-line  measure  of  the  rate  of  change  of 
normal  acceleration  would  appear  beneficial  for  reducing  bumpiness.  This  possi¬ 
bility  has  been  discussed  and  the  associated  use  of  phase-advance  filters  in 
the  feedback  path  of  the  DLC  loop  and  of  nonlinear  control  have  been  considered. 

In  addition  to  various  types  of  closed-loop  control,  this  Report  includes 
an  investigation  of  the  'open-loop*  solution.  The  latter  necessitates  the  use 
of  incidence  sensors  such  as  wind  vanes,  in  addition  to  rate  gyros  and  accelero¬ 
meters,  for  a  practical  solution. 

Comparisons  are  made  between  the  improved  ride  due  to  DLC  and  that  associ¬ 
ated  with  an  increase  of  wing-loading  for  several  aircraft.  The  objective  or 
performance  function  is  subsequently  discussed  in  detail.  The  possibility  of 
designing  a  DLC  system  that  not  only  reduces  the  rms  of  normal  acceleration  due 
to  turbulence  but  also  does  not  increase  the  g-exceedance  rate  at  low  amplitudes 
is  considered.  For  linear  closed-loop  systems  these  two  criteria  are  shown  to  be 
conflicting.  More  generally,  the  possibility  of  defining  a  'g-exceedance  level' 
curve  that  is  governed  by  both  structural  and  physiological  constraints  is 
discussed.  Some  limitations  of  the  'open-loop'  solution  are  outlined  and  the 
influence  of  rate-limiting  of  the  DLC  actuator  on  performance  and  stability  has 
been  described.  The  relative  merits  of  open-  and  closed-loop  designs  have  been 
compared. 

The  study  described  deals  with  the  rigid-body  motion  only  and  does  not 
include  effects  of  flexible  modes  of  the  aircraft.  For  a  fighter-type  aircraft, 
effects  of  structural  flexibility  would  probably  have  relatively  little  direct 
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effect  on  the  ride,  although  there  is  a  possible  structural  fatigue-load  problem 
due  to  the  decrease  in  the  tuned  gust  length  introduced  by  the  ride-smoothing 
system.  The  principles  and  techniques  involved  are,  however,  applicable  to  a 
more  complete  model  of  the  aircraft  including  flexible  modes.  The  inclusion  of 
such  modes  would  allow  possible  instabilities  due  to  coupling  between  the  control 
system  and  the  flexible  structure  to  be  monitored.  In  addition,  a  logical 
extension  planned  for  future  work  is  to  investigate  gust-load  control.  Neverthe¬ 
less,  it  is  concluded  from  the  present  study  that  whilst  possible  instabilities 
associated  with  flexible  characteristics  may  impose  constraints  on  the  degree  of 
control  in  the  case  of  combat  aircraft,  much  can  be  achieved  through  the  modifi¬ 
cation  of  the  rigid-body  modes  to  alleviate  ride  problems. 

2  DEFINITION  OF  AIRCRAFT  AND  CONTROL  SYSTEM  MODELS 

A  mathematical  model  of  a  modern  combat  aircraft  was  used  for  the  analysis. 
The  small  perturbation  equations  of  motion  are  defined  in  Appendix  A,  and  the 
aerodynamic  data  for  the  two  flight  cases  considered,  0.3  and  0.9  Mach  number, 
are  defined  in  Table  I.  Two  motivators  were  assumed  to  be  available:  elevator 
and  a  form  of  direct  lift  control  acting  via  split  flaps.  Closed-loop  feedbacks, 
pitch  rate  to  elevator  and  a  variety  of  feedbacks  to  the  DLC  motivator,  were 
considered. 

For  the  initial  investigation  simple  first-order  approximations  td  the 
elevator  and  DLC  (flap)  actuators  were  used  and  a  block  schematic  of  the  system 
is  shown  in  Fig  1.  This  simple  system  is  explored  in  depth  and  the  results 
then  applied  to  a  more  comprehensive  simulation  including  sensors,  noise  filters, 
actuators,  power  controls,  nonlinearities,  etc. 

The  mathematical  model  of  the  aircraft  contains  only  the  rigid— body  modes, 

and  does  not  include  any  flexible  modes,  although  the  implications  of  possible 

structural  coupling  are  discussed.  There  is  clearly  a  need  in  the  final  design 

of  any  active  control  system  to  include  a  fully  representative  model  of  the  air- 
11-13 

craft  including  at  least  the  first  few  flexible  modes  at  the  lower 

frequencies. 

3  ATMOSPHERIC  TURBULENCE  MODELS 

Three  different  models  of  vertical  turbulence  were  used.  They  were 
(i)  band-limited  white  noise,  (ii)  the  standard  Dryden  gust  spectrum  and  (iii)  a 
statistical  discrete  gust  model.  These  are  defined  in  Appendices  B  and  C  and  a 
non-dimensional  comparison  of  the  power  spectral  density  functions  with  the 
Von  Karman  spectrum  is  given  in  Fig  2. 

*  Subsequent  work  on  a  combat  aircraft,  however,  has  shown  that  flexibility  can 
have  an  appreciable  effect  on  the  rms  of  normal  acceleration. 
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It  can  be  shown  that  the  variance  of  a  state  output  to  white  noise  can  be 

calculated  analytically  from  an  impulse  or  initial  condition  disturbance. 

Similarly  the  variances  due  to  band-limited  white  noise  (BLWN)  can  also  be 

found  from  a  single  time  response,  by  shaping  the  initial  impulse This 

is  equivalent  to  adding  another  first  order  differential  equation  to  the  system 

model.  The  equivalent  Dryden  type  turbulence  model  can  be  defined  by  the  addi- 

1 4 

tion  of  two  first-order  differential  equations  .  Thus,  the  variance  or  rms 
of  any  output  state  to  either  BLWN  or  Dryden  type  input  can  be  calculated 
analytically  in  the  time  plane.  The  Dryden  model  was  preferred  to  the  Von  Karman 
model  simply  on  the  grounds  of  analytical  expediency,  the  Dryden  power  spectra 
being  amenable  in  particular  to  the  inverse  Fourier  transform  for  use  in  the 
state  space  domain  (see  Appendix  B) . 

3-5 

The  statistical  discrete  gust  approach  has  been  proposed  as  an  alter¬ 
native  method  for  the  assessment  of  aircraft  systems.  The  more  intense  energy 
concentrations  in  atmospheric  turbulence  are  modelled  by  equi-probable  families 
of  discrete  ramp  gusts  covering  a  wide  range  of  gradient  distance  and  intensities 
(Appendix  C) .  The  method  has  been  applied  in  this  Report  to  the  analysis  and 
assessment  of  ride-smoothing  systems  and  has  also  been  included  in  an  optimising 
procedure  as  part  of  a  digital  computer  program  (involving  minimising  the 
maximum  response  to  the  worst  gust) .  Although  slightly  unwieldy  for  digital 
computer  optimisations  the  technique  has  proved  useful  for  nonlinear  systems, 
where  the  standard  techniques  for  minimising  variances  cannot  be  applied. 

4  INITIAL  RIDE-COMFORT  OBJECTIVE  FUNCTION 

For  the  initial  phase  of  this  work  ride  improvement  is  defined  as  a 
reduction  of  the  normal  acceleration  response  of  the  aircraft  to  vertical  turbu¬ 
lence  at  specified  positions  along  the  fuselage.  This  may  be  quantified  in  two 
ways,  (i)  statistically,  by  reducing  the  rms  of  normal  acceleration  response  to 
the  standard  power  spectral  representation  of  turbulence,  and  (ii)  deterministic- 
ally,  by  reducing  the  peak  amplitude  of  the  normal  acceleration  response  to 
discrete  gusts.  The  second  method  involves  the  use  of  a  set  of  discrete  ramp 
gusts,  as  described  in  Appendix  C,  such  that  the  response  to  the  worst  gust 
can  be  investigated.  Both  methods  were  employed  and  the  results  compared. 

A  further  objective,  pertinent  to  all  systems  investigated,  was  that  all 
modes  of  the  system  should  have  damping  greater  than  a  prescribed  minimum  level. 
The  value  used  in  this  study  was  0.5  relative  damping.  As  will  be  seen  later, 
the  achievable  'ride'  improvement  is  a  function  of  this  damping  constraint. 
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Alternative  objective  functions,  other  than  merely  reducing  normal- 
acceleration  response,  are  discussed  and  analysed  in  sections  8  and  9. 

5  DESCRIPTION  OF  DESIGN  METHODS 

The  mathematical  model  used  has  been  outlined  in  section  2.  It  is  defined 
by  the  equations 

x  =  [A]x  +  [B]u  +  [C]w  (I) 

-g 

y  =  [d]  x  +  [e]  w 

L  ~g 

where  x,  u  and  w^  are  respectively  state-space,  control,  and  input  vectors 
and  y  is  the  output  vector. 

Extensive  use  has  been  made  of  an  existing  suite  of  digital  and  hybrid  (analogue/ 

digital)  computer  programs  to  design  the  feedback  control  laws  of  Fig. I.  These 
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programs  include  matrix-Riccati  optimal  control  ,  parameter  optimisation  *  , 

9  8  8  8  8  8 

frequency  responses  '  Bode  plots  ,  root  locus  ,  eigenvalues  ,  time  responses 

g 

to  deterministic  inputs  and  variances  to  statistical  inputs  .  The  parameter 
optimisation  was  the  dominant  program  used.  Given  a  form  of  control  law,  this 
program  automatically  selects  the  set  of  parameters  (feedback  coefficients,  time 
constants  etc)  which  will  minimise  a  chosen  performance  function.  The  program 
allows  for  constraints  to  be  put  both  on  parameters  and  on  cost  function  elements. 
The  variances  of  all  aircraft  states,  control  displacements  and  rates,  due  to 
the  various  turbulence  inputs,  were  calculated.  The  cost  or  performance  function 
was  a  weighted  sum  of  these  variances,  together  with  a  damping  constraint  from  all 
the  eigenvalues.  In  choosing  the  cost  function,  most  weight  was  placed  on  the 
variances  of  normal  acceleration,  pitch  rate,  and  DLC  motivator  rate. 

Another  design  method  was  to  calculate  the  normal  acceleration  responses  to 
a  family  of  discrete  gusts  (Appendix  C)  and  to  minimise  the  maximum  amplitude  due 
to  the  worst  gust. 

Further  optimisations  were  performed  in  which  the  variances  of  elevator 
and  flap  rate  were  constrained  to  within  prescribed  values. 

These  methods  were  used  to  determine  the  'optimum'  state  feedbacks  and  gain 
values  for  the  reduced  or  simple  model  (Fig  1).  This  set  the  scene  and  guide¬ 
lines  for  the  more  practical  solutions  subsequently  applied  to  specific  aircraft. 

A  block  schematic  of  the  design  procedure  is  given  in  Fig  3. 


045 


A  typical  output  of  the  parameter-optimisation  program  is  given  in  Table  2. 
The  aerodynamic  data  is  followed  by  the  initial  parameter  array  including  cost- 
function  weights  and  state  constraints.  The  variances  of  all  states  are  then 
given  for  this  parameter  set,  0 ( i )  ;  i  *  !(l)n  ,  where  n  equals  the  number  of 
parameters,  usually  the  initial  0(i)  =  0.0  .  The  variances  for  the  derivatives 
of  the  states,  eg  u,  w,  4*  9>  may  be  included.  The  system  eigenvalues 

follow,  demonstrating  the  basic  aircraft  short  period  and  phugoid  modes,  the  two 
actuator  poles,  and  finally  the  inverse  of  the  input  time  constant  (1.0  second). 

Each  of  the  N-stages  of  the  optimisation  procedure  (Table  2b),  gives  the  itera¬ 
tion  number  the  value  of  the  total  cost,  and  the  M-parameters'  values.  The 
optimisation  hill  climbing  procedure 1 finishes  either  as  a  function  of  the 
error  between  succeeding  cost  values  or  on  the  number  of  iterations.  The  end 
product  (Table  2c),  comprising  the  optimal  parameters  with  associated  variances, 
constraints  and  eigenvalues  are  then  listed.  In  this  purely  typical  example  of 
a  three  parameter  optimisation,  the  variances  of  aircraft  states  have  been 
equally  weighted  with  a  unity  factor,  and  the  state  combination  (0  -  w)  by  a 

A  I  A 

factor  of  10.0.  Here  w  =  w/V  ,  and  thus  (0  -  w)  is  a  non-dimensional  measure 
of  normal  acceleration.  No  weights  have  been  attached  to  the  variances  of  the 
derivatives  of  the  states,  ie  these  are  not  included  in  the  cost  function,  and 
there  are  no  constraints  on  any  of  the  state  variances.  If  a  constraint  had  been 
used,  associated  variances  would  have  been  included  in  the  cost  function  if 
the  current  value  was  greater  than  the  constraint.  In  this  way,  for  example, 
the  control  activity  (o’)  can  be  kept  within  prescribed  limits. 

O 

All  variances  are  in  non-dimensional  units  with  respect  to  unit  variance 

input  (w  ).  It  should  be  noted  that  each  variance  is  factored  by  the  square  of 
S 

the  weight  associated  with  it. 

For  the  'optimum'  in  Table  2c  all  modes  are  damped  to  0.5  relative  damping 
or  greater.  All  variances  have  been  reduced,  some  more  than  others,  eg  pitch 
rate  variance  has  dropped  from  0.5373  to  0.1041  and  normal  acceleration  from 
0.2931  to  0.2018.  This  relatively  small  reduction  in  normal  acceleration 
variance  is  due  to  the  choice  of  weighting  vector;  in  this  particular  case  the 
weights  on  speed  (u)  and  attitude  (0)  were,  such  as  to  put  more  emphasis  on  the 
phugoid  mode  than  necessary. 

6  RESULTS  FOR  SIMPLE  SYSTEM 

6. 1  Pitch-rate  feedback  to  elevator 

Without  the  direct-lift  motivator  operative,  and  using  only  elevator  to 
alleviate  normal  acceleration,  it  soon  became  clear  that  large  reductions 
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in  rms  normal  acceleration  at  the  pilot  station  could  be  achieved  by  feeding 
back  only  pitch-rate  to  the  elevator.  Since  the  centre  of  rotation  of  the  air¬ 
craft's  short  period  mode  is  usually  well  forward  of  the  nose,  the  pilot  will 
sense  an  acceleration  due  to  rate  of  change  of  pitch,  ie  a^  +  Iq/g  .  Both 

terms  are  reduced  by  the  pitch  autostabiliser,  due  to  the  increased  damping  of 

the  system.  Clearly,  the  less  stable  the  short  period  mode,  the  larger  the 
effect  an  elevator  loop  can'  have  on  the  ride  performance. 

Some  results  are  shown  in  Table  3  giving  percentage  reductions  in  rms 

pitch  rate  of  71.6  and  73.6,  and  in  normal  acceleration  of  17.9  and  7.4  for  the 

two  heights  of  200  and  1000  feet  respectively. 

Thus,  use  of  elevator  alone  can  produce  a  significant  element  of  ride 
improvement*.  The  system  eigenvalues  and  rms  values  of  the  other  states  in 
response  to  band-limited  white  noise  ( BLWN )  are  also  given  in  Table  3. 

A  root-locus  of  (pitch-rate  to  elevator  feedback  gain)  and  frequency- 

response  functions  (Bode  plots)  for  vertical  turbulence  input  are  given  in 
Figs  4,  5a&b.  From  the  Bode  plot  of  pitch-rate.  Fig  5a,  it  can  be  seen  why 
the  autostabiliser  reduces  the  rms  pitch  rate  so  drastically.  In  the  case  of 
the  normal  acceleration  Bode  plot,  Fig  5b,  the  main  peak  has  been  reduced  at  the 
expense  of  a  slight  increase  in  response  at  the  lower  frequencies. 

6.2  Normal  acceleration  feedback  to  direct-lift  motivator 


Rearranging  the  second  equation  of  Appendix  A 
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Thus  for  fixed  x  ,  increasing  Ba  will  always  decrease  h(a  )  (see  Fig  1). 
z  z  z 

Normal  acceleration  feedback  to  the  DLC  motivator  can  thus  be  likened  to  a 
change  in  wing-loading  (W/S),  since  t  =  m/pSV  =  W/S^pgV  .  It  should  be  noted 

*  Further  reductions  in  the  rms  of  normal  acceleration  can  be  achieved  by  feeding 
back  incidence  to  the  elevator,  but  at  the  expense  of  increases  in  the  rms  of 
pitch  rate  and  attitude. 


that  this  is  not  quite  the  same  as  reducing  the  aerodynamic  derivative  z  ,  due 

w 

to  the  influence  of  the  other  term  in  the  numerator  z  n  when  the  elevator  loop 

n 

is  closed. 


Values  for  the  two  feedback  parameters  G  and  Ba  (Fig  1)  were  obtained 

q  z 

for  a  speed  of  1000  ft/s  (M  =  0.9),  with  turbulence  inputs  applicable  to  1000  feet 

and  200  feet  height,  using  both  Dryden  and  BLWN  turbulence  models  with  cut-off 

2 

frequencies  appropriate  to  height  and  speed  . 

A  weighted  sum  of  the  variances  of  aircraft  states  was  minimised,  with  a 
relative  damping  constraint  of  0.5.  Neither  actuator  rates  or  position  were 
used  in  the  cost  function  for  this  initial  exercise. 


A  root-locus  of  Ba  with  fixed  G  (0.464)  is  given  in  Fig  6a.  Large 
z  q 

values  of  B^  can  be  employed  without  running  into  stability  problems.  The 
elevator/pitch  pole  remains  invariant  with  Ba^  ,  whereas  the  heave  pole  (real 
root)  due  to  the  DLC  loop  increases  with  B^  .  Thus  the  pitching  and  heaving 
components  of  the  aircraft  motion  have  become  de-coupled.  (An  assumption  made 
here  is  that  the  pitching  moment  due  to  the  DLC  motivator  has  been  eliminated.) 
For  the  analysis  of  this  Report  it  has  been  assumed  that  this  pitching  effect 
can  be  effectively  cancelled  by  feeding  a  signal  across  from  the  DLC  loop  to 
the  elevator  loop  with  a  gain  of  ;  ie  n  •  -  M ,/M.) 

The  pole  associated  with  the  DLC  feedback  loop  is  equivalent  to  a  rapidly 

damped  exponential  (..  The  eigenvalues  for  various  values  of  Ba^  are  given 

in  Table  4a.  A  typical  set  of  eigenvalues  with  and  without  control  are  given  in 

Table  5,  together  with  the  rms  values  of  all  states  and  motivators.  For  example, 

comparing  the  basic  and  controlled  aircraft,  the  rms  values  for  pitch-rate  (q) 

and  normal  acceleration  (a  )  have  dropped  from  2.404  and  1.127  to  0.772  and  0.254 

z 

respectively.  However,  the  large  value  of  feedback  gain  Ba^  =  -  63.5  has  also 
resulted  in  large  motivator  rate  activity:  rms  flap  deflection  of  0.122  rad,  and 
flap  rate  of  2.644  rad/s;  ie  small  displacements  but  large  rates. 

The  results  show  that,  subject  to  the  mathematical  model  assumed,  without 
any  constraints  on  the  actuator  rate  rms  the  reduction  in  normal  acceleration  rms 
is  effectively  limitless.  This  is  clearly  shown  in  ^ig  7  where  the  rms  of  normal 
acceleration  is  plotted  against  Ba^  .  The  rms  of  flap  rate  is  shown  in  Fig  8, 
also  against  Ba^  .  The  rms  of  rate  of  change  of  normal  acceleration  is  shown 
in  Fig  9.  This  is  effectively  the  only  aircraft  state  rms  that  increases  as  a 
function  of  the  feedback  gain  Bg^  .  The  significance  of  this  will  be  discussed 
later  (sections  8  and  9). 
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Figs  7,  8  and  9  also  show  the  effect  of  varying  the  actuator  time  constant 


(t2) 


where  6  = 


1  +  t2S  z  ’ 


and  from  these  it  would  appear  that  large  reductions  in  rms  aa can  be  achieved 
despite  large  actuator  lags  and  with  small  actuator  rates,  but  with  different 
values  of  the  feedback  gain  (B^),  eg  30%  reduction  in  rms  normal  acceleration 
can  apparently  be  achieved  with  time  constant  and  Ba^  values  of  0.01  and  5.0, 
and  0.3  and  12.1  respectively.  Unfortunately,  however,  with  a  more  realistic 
representation  of  the  actuation  system,  when  the  rate  limit  is  taken  to  be  a 
function  of  the  actuator  natural  frequency,  the  above  idealised  result  does  not 
hold.  This  is  described  more  fully  in  section  7.  Fig  7  also  shows  that  with  a 
fixed  gain  the  improvement  in  ride  is  a  direct  function  of  the  speed  of  response 
of  the  actuator  (t2). 


The  two  Figs  10  and  11  show  the  variation  of  o  ,  o.  ,  aj  and  o;  for 

q  az  z  6 

two  different  heights,  namely  1000  feet  and  200  feet.  Both  sets  of  results 
incorporate  a  fixed  elevator  loop  gain  of  0.464,  ie  the  zero  Ba^  values  include 
the  elevator /pitch  rate  closed  loop.  The  dotted  lines  give  the  30,  40  and  50% 
reduction  levels  for  rms  of  normal  acceleration.  It  is  interesting  to  compare 
the  results  at  the  two  heights,  the  difference  being  associated  with  the  standard 
turbulence  models,  both  Dryden  and  Von  Karman  spectrum  forms  predicting  increased 
energy  at  higher  frequencies  closer  to  the  ground. 


The  pitch  rate  to  elevator  loop  has  a  bigger  effect  on  the  rms  of  normal 
acceleration  at  the  low  altitude,  ie  a  17%  reduction  as  compared  with  an  8% 
reduction  at  1000  feet.  For  comparison  the  variations  of  rms  normal  acceleration 
for  both  heights  are  plotted  against  wing-loading  in  Fig  13.  The  effect  is  more 
marked  for  the  low  altitude  case.  However,  an  associated  difference  between  the 
two  heights  is  the  much  larger  flap  rate  (a*)  needed  to  achieve  the  same  reduc¬ 
tion  in  rms  normal  acceleration  at  the  lower  altitude.  The  two  Figs  10  and  11 
give  an  indication  of  the  actuator  characteristics  needed  to  achieve  a  signifi¬ 
cant  reduction  in  rms  normal  acceleration  (say  30%).  It  appears  that  a  minimum 
rms  of  about  10  deg/s  (actuator  rate)  is  required.  Since  the  above  graphs  are 
for  inputs  of  unity  rms  (1  m/s),  then  given  an  actuator  rate  limit  of  ±60  deg/s, 
95%  of  the  response  to  turbulence  (±3a)  would  be  within  the  actuation  rate 
limits  for  rms  input  values  up  to  2  m/s. 


Thus  it  is  clear  that  the  rate  limit  of  the  DLC  actuator  imposes  an 
important  constraint  on  the  achievable  ride  improvement.  The  size  of  the 
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feedback  loop  gain  (Bazz<j)  has  a  similar  effect.  This  gain  emphasises  the  sig¬ 
nificance  of  the  motivator  effectiveness  (Z.)  ,  since  for  large  values  of  Z  , 

o  <5 

the  equivalent  reduction  in  rms  could  be  achieved  with  lower  feedback  gain 

B3z  ,  and  hence  inferior  actuators.  Thus  the  system  performance  is  constrained 
by  a  combination  of  the  DLC  surface  effectiveness  and  the  actuation 
character i s  t i c  s . 

Fig  12  shows  the  rms  normal  acceleration  response  for  various  positions 
along  the  fuselage  for  three  cases,  basic  aircraft,  elevator  loop  only  and 
ride-smoothing  system  using  DLC.  Clearly,  with  this  aircraft,  in  order  to 
improve  the  response  at  all  stations  along  the  fuselage,  it  is  perfectly  adequate 
to  minimise  the  normal  acceleration  at  the  centre  of  gravity. 

The  effect  of  changing  the  wing  loading  of  the  basic  aircraft  is  shown  in 
Table  6  and  Fig  13.  It  can  be  seen  that  increasing  the  wing  loading  decreases 
the  rms  of  normal  acceleration  for  both  heights. 

6.3  Incidence  feedback  to  direct-lift  motivator 

Rearranging  the  second  equation  of  Appendix  A 
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Thus,  for  small  values  of  ,  incidence  feedback  is  effectively  equivalent  to 
changing  the  lift-curve  slope  z 

w 

The  eigenvalues  for  the  system  with  incidence  feedback  are  shown  in 
Table  4b  for  a  fixed  elevator/pitch-rate  loop  gain.  As  will  be  seen,  the  short- 
period  oscillatory  mode  remains  invariant  with  DLC  loop  gain.  Unfortunately, 
in  contrast  to  the  normal  acceleration  feedback,  there  is  now  a  stability  problem 
for  larger  DLC  loop  gains,  the  phugoid  mode  going  unstable  at  a  gain  value 
between  -20  and  -30.  The  root  locus  of  the  poles  is  closely  akin  to  changing  zw 
as  shown  by  the  Table  4c.  The  problem  of  phugoid  instability  is  not  serious 
since  it  can  be  solved  by  either  a  throttle/speed  control  or  a  pitch  attitude/ 
elevator  control.  Roughly  equivalent  results  to  those  obtained  with  normal 
acceleration  feedback  can  be  obtained  as  demonstrated  in  Table  5  which  shows  rms 


values  due  to  Dryden  turbulence.  A  plot  of  the  aircraft  and  control  state  rms 
values  against  is  given  in  Fig  14  for  two  values  of  DLC  actuator  time 

constant  (t^  =  0.1,  0.2).  The  trends  are  the  same  as  for  normal  acceleration 
feedback,  except  for  the  stability  boundary  at  —  27.0  . 

Although  the  theoretical  results  indicate  that  a  normal  accelerometer 
should  give  slightly  superior  'ride'  control  than  an  incidence  vane,  the  sensor 
imperfections,  noise  levels,  sensitivity  to  damage,  etc  should  all  be  considered 
in  deciding  the  choice  of  method  in  practice. 

6.4  Discrete  gust  input 

4 

This  method  involved  the  use  of  a  family  of  discrete  gust  inputs  to 
obtain  the  parameters  of  the  system  giving  the  best  performance.  The  largest 
amplitude  normal-acceleration  response  to  the  worst  gust  was  minimised.  The 
main  parameters  associated  with  this  discrete-gust  approach  are  H  ,  the  gust 
length,  y(H)  ,  the  amplitude  of  the  peak  response  and  A  ,  the  gust  length 
sensitivity,  as  outlined  in  Appendix  C.  H  and  y  are  designated  the  'tuned' 
gust  length  and  the  'tuned'  response  respectively.  Including  a  damping  constraint 
meant  that  only  isolated  ramp  gusts  need  be  used  in  this  analysis,  since  there 
were  effectively  no  significant  overswings  in  the  transient  responses. 

For  the  linear  systems  involved  the  optimum  feedback  gains  were  approxi¬ 
mately  the  same  as  for  the  minimising  variances  method.  However,  the  method, 
although  rather  unwieldy  as  a  design  technique,  proved  extremely  useful  not  only 
for  nonlinear  systems  where  eigenvalues,  etc  cannot  be  calculated,  but  as  a  means 
of  explanation  and  interface  between  different  disciplines. 

Fig  15  shows  the  maximum  amplitude  response  y(H)  of  normal  acceleration 
to  a  set  of  discrete  gusts  if  variable  length  (H) .  The  tuned  gust  length  H 
for  the  basic  aircraft  and  also  for  the  elevator-controlled  system  is  about 
70  metres.  It  may  be  seen  from  Fig  15  that  the  DLC  loop  reduces  the  peak  normal 
acceleration  response  y(H)  ,  and  that  in  addition  the  tuned  or  worst  gust 
length  H  reduces  as  the  feedback  gain  increases.  A  comparison  between  power- 
spectral-density  (PSD)  and  statistical  discrete  gust  (SDG)  methods  shows  that  a 
close  relationship  exists  between  zero  crossing  rate  ,  evaluated  using  the 
PSD  method,  and  the  quantity  1/AH  derived  from  the  SDG  method.  If  the  tuned 
gust  length  decreases  increases. 

The  'ideal'  control  system  would  be  one  that  reduced  the  amplitude  of 
response  y'H)  without  at  the  same  time  reducing  the  'tuned'  gust  length  H  . 

The  problem  of  the  changing  tuned  gust  length  is  discussed  more  fully  later  in 
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the  paper.  As  there  is  some  uncertainty  as  to  the  practical  importance  of  this 
effect,  there  is  a  need  for  further  experimental  work  to  determine  what  ride 
characteristics  the  pilot  prefers. 


7  EFFECT  OF  ACTUATOR  CHARACTERISTICS,  NATURAL  FREQUENCY,  RATE  AND 
POSITION  LIMITS 

It  has  been  shown  previously  (section  6.2)  that  the  actuator  rate  limit 
plays  an  important  part  in  determining  how  much  improvement  in  ride  can  be 
achieved.  Fig  16  demonstrates  the  effect  of  constraints  on  the  maximum  rms 

DLC  actuator  rate  (a*)  for  the  system  as  defined  in  Fig  1,  ie  normal  acceleration 

<5 

feedback  and  a  first-order  actuator.  Values  for  actuator  time  constant  are  given 
along  each  curve  of  constant  a*  .  It  should  be  noted  that  these  are  rms  values 
of  motivator  rate  per  unity  (1  m/s)  rms  amplitude  of  vertical  turbulence.  Thus 
assuming  a  Gaussian  distribution  for  the  turbulence,  (with  a  desired  physical 
limit  of  3o*),  95%  of  the  actuator  rate  would  be  within  the  linear  range. 

It  can  be  seen  that  for  a  fixed  at  ,  higher  attenuation  in  rms  normal 

o 

acceleration  is  achieved  with  more  sluggish  actuators  (larger  )  and  corres¬ 

ponding  increases  in  the  magnitude  of  the  loop  gain  Ba ^  . 

For  a  fixed  loop  gain  (B3z.) ,  further  reduction  in  rms  normal  acceleration 
can  be  achieved  by  increasing  the  allowable  rms  actuator  rate.  It  would  appear 
that  the  optimal  configuration  would  be  an  actuator  with  a  large  allowable  rms 
actuator  rate,  together  with  a  large  actuation  lag.  This  implies  large  feedback 
gains,  which  when  using  a  more  realistic  representation  of  the  actuator,  is  not 
possible  due  to  stability  constraints. 


In  order  to  elaborate  this  point,  the  actuation  system  has  been  modelled 
as  a  first  order  lag  (representing  the  power  control)  followed  by  a  second-order 
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A  root  locus  for  fixed  values  of  t„  =  0.04  second,  u>  =  50  rad,  £  =  0.7  , 

2  n 

is  given  in  Fig  6b.  Again,  the  pitching  and  heave  motions  are  decoupled,  the 
elevator-pitch  pole  remaining  invariant  with  Baz  at  approximately  -10.5  ±  18. Oi. 
However  the  interaction  between  poles  in  the  DLC  loop  causes  the  high  frequency 
mode  to  go  unstable  at  a  gain  modulus  of  just  over  50.  Similar  root  locii  occur 
for  different  values  of  u  and  £  .  The  elevator  loop  remained  fixed  to  give 
a  damping  of  0,5  for  the  pitching  mode.  A  typical  set  of  eigenvalues  for 
=  50  are  given  in  Table  4d. 
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Fig  17  gives  the  rms  of  normal  acceleration  to  unit  rms  vertical  gust 

input  plotted  against  the  natural  undamped  frequency  of  the  actuator  . 

With  this  system  the  DLC  loop  gain  is  constrained  by  stability.  The  full  line 

gives  the  0.5  relative  damping  boundary,  the  figures  in  brackets  giving  the 

relative  damping  of  the  least  damped  mode  in  the  system.  As  the  actuators 

improve  (faster  acting,  ie  higher  undamped  natural  frequency)  larger  reductions 

in  rms  normal  acceleration  can  be  achieved  in  association  with  the  larger 

gains,  Ba  .  With  increasing  gain,  the  variance  for  fixed  w  drops  until 
z  n 

the  system  approaches  the  stability  boundary  when  the  variance  increases 
rapidly. 

The  dotted  lines  in  Fig  17  show  the  30,  40  and  50%  reductions  in  rms 

normal  acceleration.  Thus,  to  achieve  a  30%  reduction  at  w  =25.0  rad  the 

n 

closed-loop  gain  would  need  to  be  about  -12.0.  However,  the  damping  of  the 
loop  would  be  reduced  to  0.34  relative  damping.  If  a  damping  constraint  of  at 
least  0.5  is  to  be  upheld  then  the  minimum  actuator  natural  frequency  would 
need  to  be  about  60.0  rad/s  (10  cps)  . 

Fig  18  shows  the  rms  values  of  actuator  rate  plotted  against  the  loop 
gain,  B3z  .  The  rms  values  are  fairly  consistent  for  different  values  of 
natural  frequency,  except  at  large  loop  gains  when  the  system  with  low  oj 

n 

tend  to  be  lightly  damped. 

It  must  be  emphasised  that  these  results  are  for  a  fixed  value  of  flap 

effectiveness,  z.  ,  ie  the  coefficient  of  lift  obtainable  from  a  deflection 
0 

of  the  DLC  motivator.  The  more  effective  the  surface,  the  lower  are  the  closed 
loop  gains  necessary  to  achieve  the  desired  improvement  in  ride,  and  hence  the 
less  actuator  rate. 

8  RE-CONSIDERATION  OF  RIDE-COMFORT  OBJECTIVE  FUNCTION. 

It  is  appropriate  here  to  introduce  the  concept  of  exceedance-level  curves. 
These  are  simply  plots  of  the  frequency  of  occurrence  N(x)  of  a  given  output 
state  exceeding  or  crossing  a  particular  level  x  . 

18  19 

It  can  be  shown  *  using  power  spectral  techniques  that  within  a 
turbulence  patch  model  led  by  a  Gaussian  distribution 


N(x) 


(6) 
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where 


Now 


and 


(7) 


N  (x)  =  -—  exp/ - (9) 

™  \  0y(H)  ) 

where  the  term  ot/AH  is  effectively  the  zero  crossing  rate.  This  latter 

approach  is  particularly  useful  for  investigating  the  response  of  nonlinear 
22 

systems 

Fig  19  shows  a  normal  acceleration  exceedance  curve  for  an  assumed 
Gaussian  patch  of  turbulence  input  of  unity  rms.  It  can  be  seen  that  has 

increased  with  increase  in  DLC  normal  acceleration  feedback,  whereas  the 
variance  represented  by  the  slope  of  the  curves  has  decreased.  also 

increases  as  the  actuator  lag  decreases.  In  contrast,  with  an  increase  of 
wing-loading,  both  rms  and  decrease  (see  also  Table  6). 
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With  the  'ride*  systems  designed  so  far,  the  rms  of  normal  acceleration 
has  been  reduced  but  the  rms  of  rate  of  change  of  normal  acceleration  has  always 
increased.  Thus  has  automatically  increased  (equation  (8)),  being  a  func¬ 

tion  of  the  ratio  of  the  two  rms  values.  However,  it  should  be  noted  from  the 
frequency  response  curve  of  Fig  5b  that  the  amplitude  of  the  normal  acceleration 
response  to  turbulence  at  high  frequencies  is  not  increased  by  the  ride-smoothing 
system  in  comparison  to  that  of  the  basic  aircraft.  Above  about  80  rad/s  (for 
this  particular  example),  the  response  amplitudes  for  both  controlled  and  basic 
aircraft  are  effectively  the  same.  Thus  the  increase  in  is  not  due  to  an 

increase  in  high  frequency  energy,  it  is  due  rather  to  an  increase  in  the  ratio 
in  equation  (8),  with  rms  normal  acceleration  (o^)  decreasing  and  at  the  same 
time  rms  rate  of  change  of  acceleration  (cu)  increasing. 

The  combination  of  reduced  rms  normal  acceleration  and  increased  have 

been  described  qualitatively  by  pilots  as  giving  a  'cobblestone  ride'.  Little 
is  known  of  the  effect  of  the  associated  increase  in  rate  of  change  of  accelera¬ 
tion  or  'jerkiness'  on  the  pilot.  However,  an  increase  in  implies  an 

increase  in  the  rate  of  occurrence  of  sign  reversals  in  normal  acceleration 
fluctuation  about  lg,  and  may  be  expected  to  have  some  adverse  effect  on  the 
pilot.  Also  any  increase  in  NQ  could  well  have  structural  implications 
associated  with  effects  on  airframe  fatigue  life. 

Previous  designs  of  'ride'  control  systems,  do  not  appear  to  have  con¬ 
sidered  these  'jerk'  effects. 

In  order  to  have  some  control  over  the  'jerkiness'  without  changing  the 

structure  of  the  feedback  system,  the  possibility  has  been  investigated  of 

including  in  the  cost  or  performance  function  to  be  minimised.  Thus  the 

cost  function  takes  the  form  +  Xa^j ,  where  X  is  an  arbitrary  constant. 

Such  a  system  would  not  only  be  required  to  reduce  the  variance  of  normal 

1  °£ 

acceleration  but  also  .  However,  since  =  - —  ,  the  cost  function 

0  0  TT  Ox 

would  effectively  be:- 

2 

°x  2 

V  =  —  +  Xax  (10) 

' 

[where  x  is  taken  to  be  normal  acceleration] . 

The  two  elements  of  the  cost  function  are  clearly  conflicting  since  mini- 
.  ,  2 

raising  o  will  decrease  the  second  term  but  increase  the  first  (w  >  1.0).  An 
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alternative  solution  would  be  to  minimise  the  variance  of  x  (normal  accelera¬ 
tion)  plus  a  heavily  constrained  ,  eg 

°x  +  C6(N0  “  N0C)  (,,) 

where  C  is  a  large  constant  and  6  is  0  or  1,  depending  upon  whether  N^  is 

less  than  or  greater  than  N  .  N  is  thus  the  prescribed  limit  N  should 

Uc  u 

be  allowed  to  achieve,  set  from  other  practical  considerations. 

On  the  basis  of  the  cost  function,  equation  (10),  however,  and  using  only 
normal  acceleration  or  incidence  feedbacks,  clearly  a  compromise  must  be  reached 
between  reduction  of  rms  normal  acceleration  and  allowing  Nq  to  increase. 

The  two  elements  of  the  cost  function  are  plotted  in  Fig  20a&b  for  different 
values  of  actuator  time  constant.  Thus  \  —  6  gives  a  30%  reduction  in  rms 
normal  acceleration  (t2=  0.04)  with  an  increase  in  from  1.5  (basic  aircraft) 

to  2.5.  Exploitation  of  this  cost  function  in  a  parameter  optimisation  procedure 
requires  additional  flexibility  (for  example  use  of  complex  filters  or  extra 
motion  feedbacks)  in  the  control  laws.  Two  possible  alternative  solutions  are 
discussed  in  section  9. 

9  ALTERNATIVE  RIDE-CONTROL  SYSTEMS 

9. 1  Shaping  of  feedback  signal 

It  was  thought  that  in  order  to  control  'jerkiness'  some  measure  of  the 
rate  of  change  of  normal  acceleration  should  be  used  in  the  feedback  loop. 

Clearly  no  sensor  measures  'jerkiness'  but  a  possible  solution  could  be:- 


This  is  equivalent  to  filtering  the  normal  acceleration  signal,  ie 


6  =  B 


B  +  B-- 
a  h 

I  +  — l -  x.S 

B  4 


(I  +  x4S)0  +  t2S) 


•'  ) 


(12) 
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where  x.  is  a  wash-out  time  constant,  and  x  represents  the  actuator  lag. 

4  £ 

Equation  (12)  may  be  written  in  the  form 

6  =  B  (1  +  X-S) /( !  +  x . S)  ( 1  +  x„S)  .  (13) 

a  3  /  4  2 

z  ' 

In  order  to  give  a  direct  comparison  with  the  normal-acceleration  feedback  only 
case  (Fig  10)  variances  for  the  relevant  states  are  given  in  Fig  21  for  fixed 
values  of  x^  =  0.01,  x^  *  0.1,  x^  =  0.05  and  varying  Baz  .  Comparing  the 
50%  reduction  in  rms  of  normal  acceleration,  for  the  non-filtered  and  filtered 
feedback  cases  the  following  rms  values  of  significant  parameters  were  obtained. 

o  0.048  (0.048) 

a 
z 

or  8.3  (7.0) 

n 

oj  16.0  (10.7) 

6 

Thus,  although  has  only  been  reduced  from  5.62  to  4.74,  the  DLC  motivator 

rate  rms  has  decreased  from  16.0  deg/s  per  unit  rms  gust  to  10.7.  Since  B3z, 

Br,  x.  (wash-out  time  constant)  and  x  (actuator  time  constant)  were  all 
h  4  z 

considered  important  parameters,  a  series  of  optimisations  were  carried  out 
with  the  following  cost  function :- 


deg/s^m/s(oJg)  . 


(14) 


where  a;  is  a  constraint  on  the  rms  of  motivator  rate,  k0  being  operative 
°C  5 

only  if  the  variance  of  actuator  rate  exceeded  the  constraint  value  and  is 

otherwise  zero,  k^  was  fixed  at  1.0,  and  varied  through  0.01,  0.02,  0.03 

for  three  optimisations. 

Table  7  gives  the  results  of  the  optimisations  with  the  values  for  the 
case  of  normal-acceleration  feedback  only  given  as  a  comparison,  ie  for  a  com¬ 
parable  reduction  in  rms  of  normal  acceleration.  From  the  table  it  can  be  seen 
that  for  the  42%  reduction  in  rms  normal  acceleration  (0.055),  the  addition  of 
a  phase-advance  network  reduces  from  4.07  to  3.27  and  rms  surface  rate 

from  12.3  deg/s  to  6.86  deg/s.  The  rms  of  surface  position  has,  however, 
increased  to  1.33  deg/m/s  rms  gust  input.  Nevertheless,  for  even  heavy 


turbulence  (3  m/s  rms)  the  surface  movement  would  still  be  less  than  12  for 

95%  of  the  time.  The  interesting  point  about  the  optimisations  is  that  the  gain 

B~  in  the  wash-out  loop  always  turned  out  to  be  opposite  in  sign  to  Ba  .  This 
n  z 

indicates  that  in  fact  a  composite  lag/lead  filter  and  not  a  lead/lag  is  the 
desired  filter  for  this  particular  cost  function  and  system. 

Optimisation  1  reduces  to  a  filter  equivalent  to  -27.4  (1  +  0. 16S)/(I  +  0.25S) 
and  an  actuator  time  constant  of  0.273  second,  thus  recommending  a  fairly 
sluggish  actuator.  This  result  may  be  seen  in  Fig  17  where  the  same  attenuation 
in  normal  acceleration  is  achieved  either  with  a  low  gain,  fast  actuator  or  a 
high  gain,  slow  actuator.  However,  a  fast  actuator  is  still  desirable  since  a 
'slow'  actuator  would  have  relatively  low  rate  limits  which  would  constrain  the 
degree  of  reduction  in  normal  acceleration  (Fig  16).  Thus  the  optimum  would 
appear  to  be  a  fast  actuator  system,  but  with  suitable  lags  in  the  feedback  loop. 

As  a  further  illustration  of  this  phenomenon.  Figs  22  and  23  show  conflict 
curves  of  var(az)  against  var(h)  and  var(a^)  against  var(6).  The  dotted  line 
shows  the  performance  for  normal-acceleration  feedback  only;  the  other  curves 
are  for  various  wash-out  time  constants  and  varying  the  B~  gain  (fixed 
actuator  time  constant  of  0.04  second).  Thus  all  the  points  to  the  left  of  the 
dashed  line  and  below  the  30%  reduction  in  rms  a ^  line  show  an  improvement 
compared  with  the  system  with  pure  normal  acceleration  feedback. 

As  stated  above,  the  studies  indicate  the  desirability  of  using  a  fast 
acting  actuator,  but  with  fairly  large  lags  in  the  loop,  and  with  high  gains. 
However,  with  a  more  realistic  simulation  of  an  actuation  system  (Fig  6b), 
stability  problems  arise  with  large  gains  (see  Fig  6b). 

Since  the  filtered  normal-acceleration  signal  appeared  to  give  beneficial 
effects  for  the  simple  actuator  system,  it  was  considered  worth  assessing  with 
a  more  realistic  actuator  model,  ie 


S  +  2_o)  S  +  to 
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As  a  baseline  for  comparison,  Fig  24  shows  a  root  locus  for  B3z(B“-  =  0.0), 

and  varying  actuator  natural  frequency  w  ^  =  0.05).  Table  8  gives  the 

variances  for  the  particular  set  of  parameters:  cj  =  120,  £  =  0.7,  t„  =  0.05, 

n  2 

Ba^  =  -12.0,  B--  =  0.0,  G  =  0.464.  With  these  parameter  values,  it  can  be  seen 

that  =  0.056,  which  is  a  40%  reduction  in  rms  normal  acceleration.  In  order 

to  see  what  improvement  could  be  achieved" for  the  same  attenuation  in  normal 

acceleration  response  and  constant  values  for  the  DLC  actuation  parameters,  a 

search  on  three  parameters  (Ba^,  B™  and  i^)  was  performed,  minimising  Nq 

with  a  constraint  of  0.056  g  on  oa  and  14.36  deg/s  on  o;  .  Table  8  gives 

z  o 

the  final  values  of  the  optimised  parameters  (Baz  =_20.5,  B~  =  12.8,  =  2.77) 

together  with  appropriate  rms  values.  As  seen  in  this  table,  addition  of  such 
filtering  enables  the  system  to  be  tuned  so  that  for  identical  rms  normal 
acceleration  there  are  reductions  in  N.  and  ot  .  Unfortunately,  however,  the 
root  locus  of  the  above  system  shows  that  a  very  small  real  root  of  the  system 
is  quite  sensitive  to  B~  .  This  root  can  become  positive,  although  it  is 
not  significant  enough  to  show  up  in  the  time  responses.  Although  this  root  can 
easily  be  stabilised  by  the  addition  of  an  attitude  feedback  to  the  elevator, 
this  has  the  effect  of  somewhat  reducing  the  damping  of  the  'pitch  rate'  mode,  so 
a  slight  re-adjustment  of  the  pitch-^rate  feedback  is  necessary.  To  illustrate 
this,  a  further  optimisation  of  the  five  parameters  ,  G0,  B^,  B-,  x^, 
again  restraining  the  reduction  in  rms  normal  acceleration  to  40%,  gives  the 
following  results: 


G  =  0.4  ,  G  =  0.53  ,  B,  =  -  20  ,  B--  =  13.6  ,  x,  =  3.2  , 

q  0  dz  h  4 

with  rms  values  as  shown  in  Table  8.  All  real  eigenvalues  are  now  well  removed 

from  the  axis  (te  no  stability  sensitivity  problem). 

A  similar  argument  for  the  minimisation  of  can  be  applied  to  the 

incidence  feedback  system.  Fig  25  shows  the  effect  on  the  state  and  control  rms 
values  of  the  incidence  system  with  the  addition  of  a  feedback  employing  rate 
of  change  of  incidence: 
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1 

(1  +  x2S) 


B  a  -  B  *  ot] 
a  a 


(15) 


It  can  be  seen  from  Fig  25  that  for  a  constant  B  (in  this  case  B  *  10), 

a  a 

increasing  the  magnitude  of  B*  results  in  a  reduction  in  ,  a*  and  o~  , 

u  U  o  h 

together  with  an  increase  in  oa^  .  From  this  figure,  and  Fig  26,  which 
illustrates  the  variation  in  performance  with  respect  to  a  family  of  discrete 


J 
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gusts,  it  can  be  seen  that  by  including  some  rate  of  change  of  angle  of  attack 
in  the  DLC  loop,  it  is  possible  to  minimise  any  increase  in  Nq  and  o*  for 
some  given  decrease  in  rms  normal  acceleration  with  respect  to  that  of  the  basic 
aircraft.  However,  since  these  rate  signals  are  obtained  by  differentiating  or 
phase  advancing  a  sensor  signal,  there  is  a  potential  problem  with  respect  to 
sensor  noise  which  could  invalidate  this  particular  solution. 

9.2  Nonlinear  control  using  dead-space 

In  this  section  the  effect  of  a  dead-space  in  the  normal-acceleration 
feedback  loop  is  considered.  Since  the  number  of  zero  level  crossings  of  normal 
acceleration  may  well  be  a  criterion  in  determining  ride  comfort,  it  was  thought 
that  attenuating  only  those  '  g'  responses  greater  than  a  particular  value  would 
improve  the  'ride'  in  this  respect.  The  discrete  gust  method  (Appendix  C) , 

was  used  to  assess  the  different  nonlinear  control  systems.  Fig  15  shows  the 
maximum  amplitude  response  of  a^yCH))  plotted  against  the  gust  length  H  for 
the  following  linear  systems:  basic  aircraft,  basic  plus  pitch  autostabiliser, 
basic  plus  pitch  autostabiliser  and  ride  control  system.  It  can  be  seen  that  H  , 
the  tuned  gust  length  at  which  maximum  response  occurs,  decreases  as  the  feed¬ 
back  gain  B3z  is  increased.  This  is  associated  with  an  increase  in  NQ  . 

The  effects  of  three  values  of  dead-space  ±0.1g,  ±0.2g  and  ±0.3g  are 
shown  in  Figs  27  and  28,  for  a  fixed  value  of  B_  =  -12.0  and  actuator  lag 

“Z 

=  0.05  second.  Fig  27  shows  the  rms  of  normal  acceleration  plotted  as  a 
function  of  the  rms  of  gust  input  (ow  ).  For  large  levels  of  input  the  system 
tends  to  the  linear  case  from  above,  but  for  ow  less  than  1  m/s,  the  ride 
control  system  has  less  effect,  especially  for  large  dead-space.  Fig  28  shows 
that  for  the  lower  levels  of  input  there  is  a  correspondingly  marked  reduction 
in  motivator  activity.  Thus  for  lower  levels  of  turbulence,  when  the  system 
is  unnecessary,  there  is  little  flap  activity,  but  for  large  levels  of  input  the 
full  flap  activity  (linear)  comes  into  play.  There  is  an  associated  reduction 
in  Nq  at  the  lower  levels  of  input.  Since  the  system  is  nonlinear  these  rms 
values  have  been  calculated  via  integration  of  a  time  response. 

9.3  Nonlinear  control  using  square  law 

In  the  previous  section,  a  nonlinear  control  law  that  separated  the 
potentially  conflicting  constituent  elements  of  the  ride  qualities  criteria  into 
those  desirable  for  low  levels  of  turbulence  intensity,  and  those  desirable  for 
the  higher  levels,  was  illustrated.  This  system,  however,  used  a  dead  zone  which 
could  render  the  results  sensitive  to  steady  or  slowly  varying  offsets  in  the 
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feedback  loop.  A  simpler  alternative  is  to  use  a  square  law  filter  *  whose 
output  xq(c)  is  related  to  an  input  x^(t)  through  the  differential  equation: 

xQ(t)  =  Ix^t)  -  xQ(t)  [  (x-^t)  -  Xq ( t ) ) x £  • 

The  overall  system  is  illustrated  in  Fig  29.  The  system's  response  characteris¬ 
tics  (measured  in  rms  quantities  in  response  to  a  150  second  BLWN  input)  for  a 
range  of  turbulence  intensities  are  show  in  Fig  30a&b,  together  with  characteris¬ 
tics  of  a  comparable  linear  system.  In  addition  Figs  31  and  32  illustrate  variation 
in  discrete  gust  response  parameters  with  turbulence  intensity.  As  seen,  for  low 
levels  of  turbulence  intensity,  the  aircraft  response  is  identical  to  that  of 
an  elevator  only  system,  with  virtually  no  DLC  activity.  As  the  intensity 
increases,  more  weighting  is  placed  on  the  normal  acceleration  response  magni¬ 
tude.  It  has  been  shown  in  earlier  sections  for  linear  ride  smoothing  systems 

that,  for  a  given  level  of  turbulence,  as  oa  is  decreased  both  a*  and 

z  6 

G£  tend  to  increase,  giving  consequently  an  increase  in  .  For  the  non¬ 

linear  system  described  above,  the  tuned  gust  length  H  will  tend  towards 
that  of  the  basic  aircraft  for  small  levels  of  turbulence  input  and  towards 
the  linear  control  case  for  large  levels  of  input.  Thus  the  ’cobblestone’ 
ride  effect  would  not  become  so  apparent  at  the  lower  levels  of  turbulence. 

Typical  time  responses  of  gust  input,  normal  acceleration  and  motivator 
rate  are  shown  in  Figs  in,  37  and  3d.  These  show  a  comparison  between  the  linear 
and  nonlinear  ride  control  systems  for  an  identical  turbulence  input.  As  can 
be  seen,  the  motivator  rate  has  been  reduced  without  much  adverse  effect  on  the 
normal  acceleration.  This  type  of  system  also  has  the  advantage  of  attenuating 
the  effects  of  sensor  noise. 

10  t)PEN’-L00P  RIDE  CONTROL  SYSTEMS 

The  principle  of  an  open-loop  ride-smoothing  system  is  to  generate  some 
estimate  of  the  vertical  component  of  turbulence,  and  to  move  the  DLC  as  a 
function  of  that  estimated  signal.  Through  this  approach,  the  handling  charac¬ 
teristics  are  unchanged,  unlike  those  of  the  closed-loop  systems  described 
previously  (although  any  handling  deficiencies  associated  with  such  closed-loop 
systems  can  be  removed  through  appropriate  filtering  of  the  pilot  stick  input). 

However,  practical  difficulties  arise  in  the  implementation  of  open-loop 
systems  as  adequate  gust  sensors  are  difficult  to  mechanise.  Wind  vanes  or 
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incidence  probes  measure  not  only  turbulence  but  also  components  of  the  aircraft 
motion,  and  it  is  necessary  to  correct  the  vane  output  with  aircraft  state 
outputs.  Furthermore,  the  integrity  of  the  incidence  vane  type  of  sensor  is  in 
doubt,  particularly  under,  say,  icing  conditions. 

An  estimate  of  the  gust  (wg  )  can  be  written  in  terms  of  three  sensor 

10  e 
outputs  ,  for  a  given  speed,  as 

w  =  w  -  I  (qV  -  h)dt  (16) 

8e  v  J 

where  is  the  output  of  the  vane  (m/s) 

q  is  the  output  of  a  pitch  rate  gyro  (rad/s) 

••  2 

h  is  the  output  of  a  vertical  accelerometer  placed  at  the  CG  (m/s  ). 

A  block  schematic  of  a  ride-control  system  using  such  open-loop  control  is 
shown  in  Fig  33.  A  detailed  analysis  of  this  system  is  given  in  Ref  10  which 
describes  how  an  open-loop  direct-lift  system  could  be  integrated  into  an 
existing  fly-by-wire  aircraft. 

Analysis  shows  that  the  open-loop  controller  becomes  slightly  more  complex 
than  the  closed-loop,  and  with  much  the  same  characteristics.  In  fact,  the 
theoretical  benefits  of  'open*  loop  become  lost  in  the  practical  implementation. 
Figs  34  and  35  show  'ride*  against  jerkiness  (section  8)  and  ’ride' 

against  surface  rate  activity  (o^),  for  varying  actuator  lags  and  loop  gain  ®Wg  • 
There  is  clearly  an  optimum  value  for  this  gain  where  the  system  is  equivalent 
to  the  closed  loop  normal  acceleration  feedback  system. 

1 1  DISCUSSION  OF  HANDLING  CHARACTERISTICS 

It  has  been  shown  that  closed-loop  direct-lift  control  systems  have  an 
effect  approximately  equivalent  to  increasing  the  wing  loading  (normal  accelera¬ 
tion  feedback)  or  decreasing  the  lift  curve  slope  (incidence  feedback).  In  both 
cases  the  dynamic  characteristics  of  the  aircraft  are  changed  (different  poles 
and  zeros).  The  short  period  pitching  and  normal  acceleration  response  to  both 
turbulence  and  pilot's  stick  input  are  changed.  In  the  absence  of  additional 
pilot-command  stick  shaping,  the  aircraft  having  a  reduced  sensitivity  to  gusts 
has  consequently  a  more  sluggish  normal  acceleration  response  to  stick  input. 

However,  if  the  stick  is  connected  to  both  the  elevator  and  DLC  motivator, 
then  by  the  use  of  stick  command  shaping  filters  practically  any  desired  transient 
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response  in  pitch,  pitch  rate,  incidence  and  normal  acceleration  can  be 
achieved.  The  effect  of  such  command  stick  shaping  is  that  the  zeros  of  the 
transfer  function  of  the  stick  to  aircraft  state  output  are  changed,  but  not  the 
poles. 

20  21 

The  theory  and  design  of  such  systems  are  discussed  in  other  papers  * 

It  is  suggested  that,  although  DLC  can  be  used  to  improve  the  ride  of  an  air¬ 
craft,  the  benefits  due  to  improved  handling  could  well  be  even  more  important, 
eg  for  air-to-ground  aiming,  weapon  delivery  etc.  Thus,  using  DLC  motivators, 
an  aircraft  designed  primarily  for  good  air-to-air  performance,  (low  wing-loading, 
large  lift-curve  slope)  could  both  achieve  good  turbulence  performance  (ride), 
at  low  altitude  and  have  handling  characteristics  optimised  for  specific  tasks. 

12  DISCUSSION  OF  RESULTS 

A  comparison  has  been  made  between  three  types  of  ride  control  system,  nor¬ 
mal  acceleration  feedback,  incidence  feedback,  and  the  so-called  open-loop  system. 
It  has  been  shown  that  approximately  the  same  reduction  in  rms  of  normal  accelera¬ 
tion  can  be  achieved  with  all  three  systems,  a  theoretical  reduction  of  at  least 
30%  being  the  target  level.  Similar  problems  arise  with  all  three  systems  and 
involve  basic  stability  and  actuator  characteristics  including  natural  frequency, 
damping  and  rate  limits,  increase  in  surface  rate  activity  and  associated  ride- 
jerkiness.  Thus,  choice  of  system  would  appear  to  depend  on  other  factors,  such 
as  complexity,  sensors,  reliability,  susceptibility  to  damage  etc.  On  this  over¬ 
all  basis  it  would  appear  that  the  use  of  accelerometers  is  slightly  superior. 


It  has  been  shown  that  ,  the  zero  level  exceedance  rate  of  normal 

acceleration,  increases  with  feedback  control,  giving  the  so-called  cobblestone 

ride.  However,  there  is  some  evidence  (from  a  survey  of  Service  pilots)  that 

pilots  prefer  a  ’hard’  ride  to  a  softer  response  with  fluctuation  of  larger 

W  /3C l 

amplitude.  This  survey  also  indicated  that  the  ratio  of  —  - - —  is  not  the 

S  /  3a 

sole  criterion  for  assessing  ride  quality  and  that  aircraft  handling  character¬ 
istics  and  trimmability  are  also  taken  into  account.  Recent  studies  indicate 
in  particular  that  'clean'  response  of  the  control  system  and  good  stick 
characteristics,  such  as  stick  force  per  g,  also  play  a  part  in  the  pilot's 
overall  assessment  of  ride  quality.  The  essential  criterion  is  whether  or  not 
the  pilot  and  crew  can  satisfactorily  perform  a  task  in  turbulence. 

The  paper  also  indicates  that  a  good  pitch  autostabiliser  is  always  bene¬ 
ficial  and  in  some  cases  can  reduce  the  rms  of  normal  acceleration  by  up  to  20%. 
The  effect  of  height  variations  on  the  frequency  content  and  amplitude  of  the 
turbulence  is  also  relevant  to  the  design. 
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Most  ride  systems  designed  in  the  past,  have  used  existing  aerodynamic 
surfaces  to  generate  the  direct  lift.  Clearly,  this  is  unsatisfactory,  since 
they  were  not  designed  for  this  specific  purpose,  and  more  research  is  required 
with  the  objective  of  providing  more  effective  direct-lift  devices  preferably 
with  little  or  no  pitching  moment. 

The  effect  of  ride-control  systems  on  the  handling  characteristics  is  also 
very  important  and  is  the  subject  of  another  paper.  Suffice  to  say  here,  that 
advantage  can  be  taken  of  the  extra  motivator  to  achieve  any  desired  responses 
to  pilot's  stick  in  both  pitch  and  heave. 

One  major  limitation  of  the  study  presented  in  this  paper  has  been  the 
omission  of  consideration  of  the  structural  loading  changes  and  possible 
instabilities  due  to  coupling  between  the  control  system  and  the  flexible  air¬ 
craft  structure.  There  is  thus  a  need  for  continuing  research  into  such  effects. 

The  use  of  approximate  mathematical  models  in  the  design  process  has 
proved  very  useful,  particularly  in  the  physical  understanding  of  the  different 
control  systems.  However,  there  are  dangers  in  using  simplified  models,  especi¬ 
ally  with  respect  to  high  feedback  gains  and  general  stability,  as  there  may 
exist  limitations  not  necessarily  indicated  by  the  simplified  model.  Also,  the 
practical  limitations  of  the  actuation  system  should  be  included  in  the  final 
assessment. 

The  optimisation  of  a  ride-smoothing  system  using  discrete-gust  theory 
(Appendix  C) ,  ie  minimising  the  maximum  response  to  the  worst  gust,  gave  similar 
system  parameters  as  optimisation  using  the  more  usual  power-spectral  methods 
(Appendix  B) .  Although  the  discrete-gust  method  is  rather  lengthy  and  unwieldy, 
the  present  study  indicated  that  the  technique  can  be  used  for  designing  nonlinear 
control  systems  and  in  fact  the  theory  is  currently  being  extended  for  this 
purpose. 

In  order  to  give  some  time  responses  of  state  parameters  the  analogue  model 
of  the  aircraft  plus  control  system  was  excited  by  a  record  of  actual  turbulence, 
the  results  of  which,  for  pitch  rate  and  normal  acceleration,  are  shown  in 
Figs  36-38  for  the  uncontrolled  and  controlled  aircraft  (normal  acceleration 
feedback).  The  two  normal  acceleration  responses  were  subsequently  used  to  drive 
a  motion/vibration  simulator  for  pilot  assessment.  This  latter  experiment  in 
the  human  factors  area  is  continuing. 
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13  FINAL  COMMENTS 

It  has  been  shown  that  reasonable  reductions  in  rms  normal  acceleration 
(*30%)  are  potentially  achievable  using  a  properly  designed  ride-smoothing  system 
with  direct-lift  control  (DLC) .  The  recommended  system  uses  accelerometers  to 
measure  normal  acceleration  as  the  input  signal  to  the  DLC  loop.  This  results 
in  a  system  with  comparatively  fast  motivator  surface  rates  but  relatively  small 
displacements.  Consequently,  there  is  a  need  for  fast  (high  natural  frequency) 
actuators  with  wide  rate  limits.  Although  theory  suggests  that  the  same 
rec  iction  in  rms  of  normal  acceleration  can  be  achieved  with  either  low  loop 
gains  and  fast  actuators  or  high  loop  gains  and  slow  actuators,  stability  prob¬ 
lems  can  occur  with  the  higher  loop  gain  systems.  Consequently,  a  compromise  has 
to  be  reached  between  stability,  reduction  in  rms  of  normal  acceleration  and 
surface  rate  activity.  The  optimum  system  appears,  paradoxically,  to  require 
fast  actuators  but  with  some  lag  inserted  into  the  feedback  loop. 

A  consequence  of  controlling  the  normal  acceleration  through  the  use  of 
DLC  is  that  the  rms  rate-of-change  of  normal  acceleration  increases.  This  leads 
to  an  associated  increase  in  the  'zero  crossing'  parameter  N^  (which  reflects 
the  overall  rate  of  sign  reversal  in  the  'g'  fluctuations).  This  'jerkiness'  or 
increase  in  NQ  could  well  have  a  significant  effect  on  the  quality  of  the  ride, 
and  although  there  is  some  evidence  that  pilots  prefer  this  'hard'  or  'cobble¬ 
stone'  ride,  further  human  factors  work  is  required  to  fully  evaluate  the 
implications.  This  paper  suggests  that  Nq  can  be  constrained  by  appropriate 
design  of  the  control  system.  In  particular,  work  has  indicated  that  nonlinear 
control  systems  have  promising  potential  for  this  purpose. 

Pilot  surveys  have  indicated  that  ride  quality  associated  with  the  basic 

.  .  W  / 

airframe  is  not  merely  a  function  of  the  parameter  •=■  / Cl  ,  and  that  handling 

o  /  c t 

qualities,  including  the  ability  to  control  the  aircraft  response  in  pitch, 
are  equally  important. 

The  potential  benefits  of  DLC  are  twofold:  improvement  of  the  ride  and 
also,  using  the  extra  motivator,  improvement  of  the  handling  qualities,  particu¬ 
larly  for  target  acquisition  and  tracking.  Proposed  future  work  includes  the 
assessment  of  control  systems  employing  aerodynamic  surfaces  (and/or  deflected 
< hrust)  designed  with  this  joint  objective  in  mind. 
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Appendix  B 

SPECTRAL  TURBULENCE  MODELS 


B. I  Band  limited  white  noise  (BLWN) 


White 

noise 

generator 

Linear 

BLWN 

filter 

ie  BLWN  can  be  obtained  by  passing  white  noise  through  a  linear  filter,  or  if 
variances  of  the  output  signal  are  required:- 
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ie  variances  o£  state  outputs  can  be  obtained  by  using  an  impulse  or  putting  an 
initial  condition  on  a  linear  filter.  This  can  be  expressed  as  a  first  order 
differential  equation  by:- 

x  -  "  7  x  (B-l) 


where  t  determines  the  cut-off  frequency  of  the  power  spectral  function  and 
the  initial  condition  is  of  magnitude  /2/t  to  give  unity  variance  of  the 
input . 

B.2  Dryden  spectrum 

Dryden  gives  the  power  spectral  function  for  vertical  turbulence  as: 


$ 

w  w 


g  g 


(«) 


TT 


where  L  =  gust  length 

V  =  aircraft  velocity. 

By  using  the  inverse  Fourier  transform  equation  (B-2)  gives :- 


(B-2) 


yw  (t) 
g 


/3  +  (1  -  /3)  — 

Li 


(B-3) 


In  a  similar  way  to  BLWN, 
differential  equations 


y»g<t) 


can  be  expressed  as  the  solution  to  two 


the  initial  conditions  being:- 
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y.(0)  -  y  (0) 

1  'w 

g 

y,(0)  =  y  (0) 

2  w 

g 

These  equations  can  be  used  as  an  extension  of  the  system  [A]  matrix  in  order  to 
generate  the  system  variances  to  Dryden  type  turbulence. 

B. 3  Von  Karman  spectrum 

The  power  spectral  function  for  vertical  turbulence  is  given  as:- 
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xx 
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(B-6) 


Unfortunately  this  equation  is  not  easily  amenable  to  inverse  Fourier  transform 
so  cannot  be  used  in  a  comparable  way  to  both  BLWN  and  Dryden  for  direct  calcu¬ 
lation  of  variances.  Comparisons  of  the  three  power  spectral  functions  are 
given  iq  Fig  2. 
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Appendix  C 

STATISTICAL  DISCRETE  GUST  RESPONSE  ANALYSIS 
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The  following  description  of  a  family  of  discrete  ramp  gusts,  and  the 
definition  of  aircraft  performance  in  terms  of  the  response  to  these  gusts,  are 
based  on  Ref  5, 

The  family  of  discrete  gust  is  illustrated  in  Fig  39,  and  is  defined  through 
equations  (B-l)  and  (B-2) : 


w 

g 


h  <  0 

0  ,<  h  s<  H 

h  >  H  , 


(C-l) 


uoH 


'uJLT 
0  w 


» 


Ha  l 

w  I 

[  (C— 2 ) 

h>,i  I 

w  I 


where  h  is  the  penetration  distance  through  the  gust, 

H  is  the  gust  length  or  gradient  distance, 

Uq  is  a  gust  intensity  parameter 

and  L  is  the  turbulence  scale  length, 
w 

For  well-damped  systems,  such  as  those  investigated  in  the  present  study, 
the  normal-acceleration  response  to  the  above  family  of  ramp  gusts  is  generally 
such  that  there  is  a  one-to-one  correspondence  between  significant  response  peaks 
and  individual  gusts'*  (for  less  well-damped  systems,  gust  patterns  involving 
interacting  ramp  components  need  to  be  considered). 


The  responses  may  be  assessed  in  terms  of  the  scaled  magnitude 
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of  the  amplitude  a  (H,u„)  of  the  largest  peak  in  the  normal-acceleration 
z  u 

response  to  an  isolated  ramp  gust  of  length  H  as  defined  by  equation  (C-2) . 

In  the  case  of  a  linear  system,  y  ^  (H,Uq)  becomes  independent  of  u^  and  we 
have 

Y  (H)  =  y  (H,u  )  .  (C-4) 

a  a  u 

z  z 


Examples  of  the  functions  y3z(H)  for  a  linear  system  and  Ya^Ol.u^)  for 
a  nonlinear  system  are  illustrated  in  Figs  15  and  31  respectively. 


By  varying  the  gust  length  H  it  is  possible  to  find  a  tuning  condition 
corresponding  to  a  peak  in  the  curves  Yaz(H)  or  Ya^^.u^)  .  We  denote  these 
peak  values  by 


Y  (H) 
a 


Ya  (U0} 
z 


Ya  (H(u0),u0} 
z 


and  the  associated  tuned  gust  lengths  by  H  or  H(Uq)  *  In  addition,  a  gust 
length  sensitivity  A  ,  which  is  a  measure  of  the  'breadth*  of  the  peak  in  a 
y(H)  curve,  may  be  defined'*.  For  most  practical  applications  the  variations 
in  A  are  relatively  insignificant  and  for  the  purposes  of  the  present  work  a 
nominal  value  of  A  =  0.2  may  be  assumed. 


The  discrete-gust  response  assessment  is  then  based  on  the  performance 
measures 


Y  >H  (linear  system) 

a 

z 


or 


Y  (u0) ,  H(u,,)  (nonlinear  system), 
a  u  u 

z 


In  statistical  terms,  the  rate  of  occurrence,  per  unit  distance  travelled, 
of  response  peaks  greater  than  magnitude  y  is  given  by 


N 

y 


a 

AH 


(05) 


where  a  and  8  are  turbulence-dependent  parameters.  In  Ref  5  the  values 


045 


Appendix  C 


35 


a  =  0.38 


0.075  (ft/s,  units) 

0.10a  (m/s  units) 


proposed  for  5  mile  patch  lengths  at  low  altitude,  where  a  is  a  turbulence 
reference  intensity.  Thus  a  is  assumed  to  be  a  constant  and  8  is  a  constant 
for  each  turbulence  patch  but  varies  from  patch  to  patch  in  proportion  to  5  . 

In  the  case  of  a  nonlinear  system  it  is  convenient  to  take  u^  as  an 
independent  variable.  From  equation  (C-3)  it  then  follows  using  the  tuning 
condition  that  the  corresponding  value  of  y  is  given  by 


y 


una 
0  z 


<v 


(C-6) 


Also,  equation  (C-5)  may  be  written  in  the  form 


N 


y 


a 

yH(uq) 


(C-7) 


From  equations  (C-6) 

N  and  y  with  u„ 
y  0 

By  letting 


and  (C-7)  we  thus  obtain  a  functional  relationship  between 
as  a  parametric  variable. 

tend  to  zero  we  obtain  a  quantity 


N 


0 


(C-8) 


Analogous  to  the  'zero-crossing  frequency'  (equation  (7)  in  main  text)  of 

the  power-spectral  method. 
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Table  1 

AERODYNAMIC  DATA  FOR  THE  TWO  DATUM  CASES 


Flight  case 

M  =  0.3 

1 

Flight  case  2 

M  =  0.9 

V 

103.1 

m/ s 

307.7  m/s 

X’ 

0.055 

0.040 

U 

X’ 

0.143 

0.002 

W 

k 

0.104 

0.032 

X6 

0.0 

0.0 

Z* 

0.208 

0.045 

U 

z' 

0.674 

1.696 

W 

ZQ 

0.012 

0.002 

Z* 

0.137 

0.246 

n 

Z6 

0.02 

0.082 

M- 

0.0 

0.281 

u 

M’ 

1.882 

21.066 

w 

Ml 

0.141 

0.670 

w 

M' 

0.539 

1.27 

q 

M' 

6.381 

42.134 

n 

M6 

0.50 

3.652 

Mass 

10000  kg 

Wing  area 

24.06  m2 

Span 

8.49  m 

W/S 

4078  N  nf2 
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Table  2a 

A  TYPICAL  PARAMETER  OPTIMISATION  OUTPUT 
INITIAL  PARAMETER  ARRAY 


Optimise 

G  >  B 

q  9 

a 

Flight  case  1 

X*  = 
u 

0.055 

,  z* 

w 

0.674  , 

M'  =  1.882 

w 

xf  * 

w 

0.143 

2  ’  = 

’  Q 

0.012  , 

m:  =  o.ui 

w 

xe  = 

0.104 

z '  = 

’  n 

0.137  , 

M’  =  0.539 

q 

X6  = 

0.0 

?*  s 

’ 

0.02  , 

M.  =  6.381 

n 

Z'  = 
u 

0.208 

,  M’ 

0.00  , 

Mj  =  0.50 

State 

Weight 

Variance* 

Constraint 

u 

1.0 

0.2005 

0.0 

w 

1.0 

0.9178 

0.0 

q 

1.0 

0.5373 

0.0 

6 

1.0 

0.8829 

0.0 

n 

1.0 

0.00 

0.0 

6 

1.0 

0.00 

0.0 

w 

g 

1.0 

1.00 

0.0 

0  -  w 

10.0 

29.3148 

0.0 

System  poles 

(eigenvalues) 

Damping 

Frequency 

Relative  damping 

-0.6846 

-0.6846 

1.3145 

-1.3145 

0.4619 

0.4619 

|  short  period 

-0.0191 

-0.0191 

0.1348 

-0.1348 

0.1403 

0.1403 

|  phugoid 

-20.0 

0.0 

1.0 

elevator  actuator 

-25.0 

0.0 

1.0 

DLC  actuator 

-1.0 

0.0 

1.0 

turbulence  input 

*  Non-dimensional  units. 


Table  2b 


TYPICAL  SEARCH  PROCEDURE 
Parameters 


Cost  value 

G 

G. 

B 

a 

0 

a 

0 

7.420777F 

01 

2.000000 

0.000000 

0.000000 

0 

2.313220E 

01 

1 .000000 

1.732051 

0.000000 

0 

2.426833E 

01 

1  .  oooooo 

0.577350 

1 .632993 

0 

3.285332E 

01 

0. oooooo 

0.000000 

0.000000 

SMALLEST  COST= 

23.132201 

1 

2.446655F 

01 

2.666667 

1 .539601 

1 .088662 

2 

2 . 397088F 

01 

2.000000 

1  .  154701 

0.816497 

SMALLEST  COSTs 

23.1 32201 

3 

2. 475312E 

01 

2.333333 

1 .347151 

-1 .088662 

4 

2.364007E 

01 

1  .  333333 

0.769800 

0.952579 

SMALLEST  COSTs 

23.132201 

5 

2 . 398765  F 

01 

0.888889 

2.437701 

1.179384 

ft 

2. 365747E 

01 

1 .166667 

1.828276 

0.884538 

SMALLFST  COSTs 

23.132201 

7 

2.269598F 

01 

0.333333 

1.732051 

0.408248 

8 

1 . OOOOOOE 

05 

-0.500000 

2.020726 

0.204124 

smallest  cost= 

22.695985 

9 

2 . 234427F 

01 

0.611111 

0.994325 

0.022680 

10 

2. 186945E 

01 

0.333333 

0.577350 

-0.408248 

SMALLEST  COSTs 

21 .869454 

11 

1 .OOOOOOE 

05 

-0.22222? 

1  .924501 

-0.952579 

12 

2.294269E 

01 

0 .9  44444 

1.058475 

0.476290 

SMALLFST  COSTs 

21.869454 

13 

2.747502E 

01 

0.074074 

0.513200 

0.317526 

SMALLEST  COSTs 

21 .869454 

14 

1 .OOOOOOE 

05 

-0.450617 

0.823259 

-0.264605 

15 

2.238345E 

01 

0.595679 

0.999671 

0.291066 

SMALLFST  COSTs 

21 .869454 

1ft 

1 . OOOOOOE 

05 

0.335391 

-0.338570 

-0.274686 

17 

2.221826F 

01 

0.333848 

1.21439ft 

0.237515 

SMALLFST  COSTs 

21 .869454 

18 

2.275418F 

01 

0.767833 

1 . 347745 

-0.237305 

19 

2.177888E 

01 

0.247514 

0.721836 

0.178819 

SMALLFST  COSTs 

21 .778879 

20 

2.458962F 

01 

0.014118 

0.676050 

-0.285676 

21 

2.708664E 

01 

0.450289 

0,918766 

0.146880 

SMALLEST  COSTs 

21 .778879 

22 

2. 220949E 

01 

0.353576 

0.264239 

-0.292548 

23 

2.181530E 

01 

0.348644 

0.501778 

-0.160032 

SMALLEST  COST= 

21 .778879 

24 

2.263796E 

01 

0.169372 

0.281877 

-0.406522 

25 

2. 1881 08E 

01 

0.380059 

0.759544 

0.008530 

SMALLFST  COST= 

21 .778879 

26 

2 . 1 8869  5  F 

01 

0.239601 

0.441100 

-0.268171 

27 

2.180596E 

01 

0. 344945 

0.679933 

-0.060645 

SMALLEST  COSTS 

21 .778879 

28 

2.182776E 

01 

0.294068 

0.691681 

0.380343 

29 

2 . 1 771 69  E 

01 

0.303885 

0.663099 

0.183195 

SMALLEST  COSTS 

21.771687 

30 

2.192830E 

01 

0.248918 

0.874800 

0.360944 

31 

21.815296 

21 

.805961  21 

.778879  21 

.771687 
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Table  2c 


FINAL  PARAMETER  ARRAY 


State 

Weight 

u 

1.0 

w 

1.0 

q 

1.0 

e 

1.0 

n 

1 .0 

6 

1.0 

w 

1.0 

8 

0  -  w 

10.0 

System  poles 

(optimum) 

Damping 

Frequency 

-25.0004 

0.0 

-1.4422 

2.2193 

-1.4422 

-2.2193 

-0.0209 

0.0 

-0.4526 

0,0 

-18.0491 

0.0 

-1.0 

0.0 

Variance 

Constraint 

0.0073 

0.0 

0.3971 

0.0 

0.1041 

0.0 

0.0397 

0.0 

0.0267 

0.0 

0.0149 

0.0 

1.0 

0.0 

20.1819 


Relative  damping 

1.0 

0.5449 

0.5449 

1.0 

1.0 

1.0 

1.0 


Optimal  parameter  values 


G 

q 

0.3034 


0 

0.6631 


B 

a 

0.1832 


in 

<r 

o 


i 


TR  7904b 


Table  4 


Poles 


Poles 


-0.032 

-0.02  ±  0.001 2i 
-0.02  +  0.01 4i 
-0.02  ±  0.023i 
-0.02  ±  0.034i 
+0.016  ±  0.0951 


=  0.464 


-10.6  +  18.01  -0.0087 

»  M 


x.  »  0.05  ,  G  =  0.464 
'  9 


Poles 


-2.50 

-0.032 

-0.0087 

-2.30 

-0.029 

-0.0120 

-1.82 

-0.02 

±0.01 I7i 

-1.34 

-0.02 

±0.0220i 

-0.86 

-0.019 

±0.0356i 

-0.11 

+0.09 

±0.087i 

10.6  ±  18.01 


t2  «  0.04  , 


-0.038 


-0.032 


Poles 


=  0.464 


-0.006 

-0.0058 

-0.0047 

-0.0028 

-10.5  ±  18.01 

tl 

II 

II 

-21.7  ±39.7i 
-18.3  ±  42. 8i 
-13.6  ±  47. 9i 
-3.88  ±  60. 7i 

-0.0087 

-10.6  ±  18. Oi 

-35.0  +  35. 7i 

Table  5 


EIGENVALUES  AND  RMS  VALUES  FOR  CONTROLLED  AND  UNCONTROLLED  AIRCRAFT 


BASIC  AIRCRAFT  M  =  0.9  HEIGHT  1000  FT 


-1.80 
-1 .80 

-0.022 

-0.022 


-1.0 

-0.997 


Frequency  1  Relative  dampin 


0.248 

0.248 


short  period 
phugoid 

elevator  actuator 
DLC  actuator 

Dryden  input 


Rms  value*  (Dryden  input) 


State 

Basic  a/c 

Normal  acceleration 
feedback 

Incidence  feedback 

u 

0.094 

0.0837 

0.109 

h 

0.356 

0.162 

0.159 

q 

2.404 

0.772  (68%) 

0.774  (68%) 

e 

0.935 

0.526 

0.618 

n 

0.0 

0.342 

0.349 

6 

0.0 

0.122 

0.141 

u 

0.030 

0.016 

0.019 

h 

1.127 

0.254  (78%) 

0.412  (64%) 

• 

12.827 

7.621 

7.290 

6 

2.404 

0.772 

0.774 

t 

0.0 

2.422 

2.348 

5 

0.0 

2.644 

1.135 

w 

1.0 

g 

Freauenc 


-158.21 
-10.468 
-10.468 
-0.022 
-0.022 
-1.213 
-1 .00 
-0.997 


0.0 

18.132 

-18.132 

0.025 

-0.025 

0.0 

0.0 

0.0 


*  Rms  values  are  in  non-dimensional  units  for  comparison  only 


Table  6 

EFFECT  OF  WING-LOADING  ON  KMS  NORMAL  ACCELERATION  AND  N 


Basic  aircraft  varying  wing  loading, 
(a)  I  second  time  constant  (1000  ft),  BLWN 


%  change 
with 
respect 
to  basic 
aircraft 

a 

a 

z 

'g's/l  m/s  gust 

-20% 

0.112 

0% 

0.094 

+20% 

0.082 

+40% 

0.072 

+60% 

0.065 

+80% 

0.059 

(b)  0.2  second  time  constant  (200  ft),  BLWN 


-20% 

0.2 

2.78 

0% 

0.167 

2.633 

+20% 

0.144 

2.539 

+40% 

0.126 

2.467 

+60% 

0.113 

2.415 

+80% 

0.102 

2.374 
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-10.53  ±  1 7.95i  -9.98  ±  16.261 

-5.22  -6.628 

-1.366  -0.042 

-0.0387  -2.295 

+0.016  -0.229 
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LIST  OF  SYMBOLS 

normal  acceleration  in  g  s 

normal  acceleration  feedback  gain  (DLC  loop) 

incidence  feedback  gain  (DLC  loop) 

rate  of  change  of  incidence  feedback  gain  (DLC  loop) 

incidence  feedback  gain,  open  loop  system 

band  limited  white  noise 

gradient  distance  of  discrete  gust 

direct-lift  control 


feedback  gains  of  elevator  loop 


'open'-loop  cross-feed  gain  to  elevator 
discrete  gust  length 
'tuned'  gust  length 
0  -  w  =  0  -  a 
0  -  w  =  q  -  4 

arbitrary  constant 

gust  length  (Dryden  spectrum) /scale  length  of  discrete  gust 

Mach  number 

(V  — 

Nq/2J  e  ,  exceedance  frequency 

zero-level  exceedance  frequency 
rate  of  change  of  pitch  attitude 
root-mean  square 
Laplace  operator 
m/pSV 

aircraft  velocity  (m/s) 
aircraft  weight 
estimate  of  gust  velocity 
output  of  wind  vane 
wing  loading 


a 

& 

Y  (H) 
Y(H) 
6 
<t> 


angle  of  incidence 

rate  of  change  of  incidence;  da/dt 

amplitude  of  first  overshoot  of  response 

'tuned'  amplitude 

DLC  motivator  deflection 

power  spectral  density  function 


o 

O' 


XX 


LIST  OF  SYMBOLS  (concluded) 

a.  rms  of  state  '  i' 

1 

Tj  time  constant  of  elevator  actuator 

■*2  time  constant  of  DLC  motivator  actuator 

t ^  time  constant  of  shaping  filter 

time  constant  of  shaping  filter 

x  wash-out  time  constant 

6  pitch  attitude 

X  weighting  parameter  in  cost  function  /discrete  gust  sensitivity  parameter 

n  elevator  angle 

C  damping  ratio 

wn  undamped  natural  frequency 

Aerodynamic  symbols  and  derivatives 
u  =  u/V  forward  speed 

A 

w  =  w/V  vertical  speed 
q  pitch  rate 

0  pitch  attitude 

A 

fi  =  0  -  w  non-dimensional  rate  of  change  of  height 
n  elevator  angle 

6  flap  (DLC)  angle 


X'  =  -  x  /t 
u  u 

X'  =  -  x  /t 
w  w 

1  A 

k  =  CL/C 

Z’  =  -  z  /t 
u  u 

Z  =  -  z  /t 
w  w 

z;  -  -  v£ 
n  -  - 

t  *  m/pSV 


■  -  •i,»/y! 

t  -  -  “imu/iB£2 

■  -  vyv 

M'  =  -  m  /it 
q  q  B 

K  ’  ‘  Uimr,nBl2 

’  -  >'lm«/iB£2 


Pj  =  m/pS^t 
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Frequency  response  of  normal  acceleration  to  gust  input 
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with  Ba, 
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Fig  8  Effect  of  feedback  gain  and  actuator  time  constant  on  rms  of  OLC 
motivator  rate 
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Fig  10  Effect  of  feedback  gain  (Ba  j  on  rms  values.  1000  ft,  7*2  =  0.04 
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Fig  14  Effect  of  incidence  feedback  on  rms  values 


Fig  15  Normal  acceleration  maximum  amplitude  response  to  discrete  gusts.  Altitude  200 
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Fig  37  Response  of  linear  RCS  (B 
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